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Abstract 

We have performed the astrometry of H2O masers in the Galactic star- forming region Onsala 2 North 
(ON2N) with the VLBI Exploration of Radio Astrometry. We obtained the trigonometric parallax of 
0.261 ±0.009 mas. corresponding to the heliocentric distance of 3.83 ±0.13 kpc. ON2N is expected to be 
on the solar circle, because its radial velocity with respect to the Local Standard of Rest (LSR) is nearly 
zero. Using present parallax and proper motions of the masers, the Galactocentric distance of the Sun 
and the Galactic rotation velocity at the Sun are found to be i?o = 7.80 ±0.26 kpc and 0o = 213 ±5 km 
s -1 , respectively. The ratio of Galactic constants, namely the angular rotation velocity of the LSR can be 
determined more precisely, and is found to be Oo = ©o/-Ro = 27.3 ±0.8 km s _1 kpc -1 , which is consistent 
with the recent estimations but different from 25.9 km s _1 kpc -1 derived from the recommended values of 
0o and Rq by the International Astronomical Union (1985). 
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1. Introduction 

Very Long Baseline Interferometry (VLBI) astrometry 
is an important method to measure the structure of the 
Milky Way Galaxy (MWG). By measuring the accurate 
position of the source and its time variation, the source 
distance and proper motion can be determined directly. 
VLBI astrometry at 10 /ias accuracy of the Galactic H2O 
and CH3OH maser sources with the VLBI Exploration 
Radio Astrometry (VERA) and the Very Long Baseline 
Array (VLBA) can determine accurate distances at kpc- 
scalc with the errors less than 10% (see for example 
Hachisuka et al. 2006; Xu et al. 2006; Honma et al. 2007). 

The galactocentric distance of the Sun, Rq, and the 
Galactic circular rotation velocity at the Sun, 0o, are 
two fundamental parameters to study the structure of 
the MWG and they are called here as the Galactic con- 
stants. The rotation curve of the MWG and all kine- 
matic distances of the sources in the MWG are derived 
from these parameters. Since 1985, the International 
Astronomical Union (IAU) has recommended to give the 
values of Rq = 8.5 kpc and 0o = 220 km s — 1 . However, 
recent studies report the values different from them (e.g. 
Miyamoto & Zhu 1998; Reid et al. 2009a). 

* Present address: Graduate School of Science and Engineering, 
Kagoshima University, 1-21-35 Korimoto, Kagoshima, 
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Estimation of the Galactic constants is, however, af- 
fected by several independent assumptions; the peculiar 
motion of the source, systematic non-circular motions of 
both the source and the Local Standard of Rest (LSR) due 
to the spiral arm potential, the non-axisymmetric poten- 
tial of the MWG, the warping motion of the galactic disk, 
and, furthermore, the motion of the Sun with respect to 
the LSR. In this paper to simplify the situation we assume 
that the source moves with the perfect circular rotation 
on the disk. 

The circle with the raduis of the galactocentric distance 
i?o of the Sun is called as the solar circle. All the sources 
on the circle perform the circular motion with the circu- 
lar velocity Oo of the LSR, provided that the non-circular 
motion of a source in the galactic disk is negligible. Due to 
the symmetric geometry, the radial velocity of a source on 
the circle is observed to be zero with respect to the LSR, 
and the proper motion of the source depends only on the 
Galactic rotation velocity Go of the LSR. Therefore, we 
can derive 0o from the measured proper motion of the 
source on the circle. We can also derive Rq from the he- 
liocentric distance of the source, since the source, the Sun, 
and the Galactic center make an isosceles triangle (Figure 
1). Thus, we can find directly the value of the angular 
velocity fio of the LSR from 8o/i?o- Traditionally, this 
value has been derived from the Oort constants A and B 
on the basis of the kinematic analysis of stars in the solar 
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Fig. 1. The method to estimate Ro and ©o from the distance 
and the proper motion of the source on the solar circle. 

neighbourhood (see Miyamoto & Zhu 1998). The value of 
this ratio is a constraint to estimate one of the Galactic 
constants from the other. Although the IAU gives recom- 
mended values of the Galactic constants, at least one of 
them should be revised, if the ratio is inconsistent to the 
observed value. 

Onsala 2 North (ON2N) is a massive star-forming 
region located at the Galactic coordinates of (I, b) = 
(75.°78,0.°34). Its radial velocity i>lsr with respect to the 
LSR is observed to be Oil km s _1 in the NH 3 and CS lines 
(Olmi & Cesaroni 1999; Codella et al. 2010). Lekht et al. 
(2006) detected the H 2 masers of ON2N at the radial 
velocity range from —12 to 9 km s _1 , with their peak flux 
densities of 10 2 -10 3 Jy. These H2O masers are associated 
with a 7 mm radio continuum source and an NH3 core as 
well, which are located at approximately 2" south from 
the ultracompact Hll region, G75.78+0.34 (Carral et al. 
1997; Codella et al. 2010). Thus, ON2N is considered to 
be one of the sources on the solar circle. 

We made astrometric observations of the H2O masers 
in ON2N with VERA. Based on these proper motions 
and trigonometric parallax measurement, we estimate the 
Galactic constants Qq, Rq, and fio- 

In this paper we use the "LSR velocity" or i>lsr as the 
radial velocity with respect to the frame moving toward 
ai900 = 18 h , <5i9oo = +30° with —20 km s _1 namely, a pro- 
visional solar motion with respect to the LSR to be deter- 
mined after the traditional definition in radio astronomy 
since 1960s (see Kerr & Lyndcn-Bcll 1986). 

2. Observations and Data Reductions 

We observed H2O masers at the rest frequency of 
22.235080 GHz in the star-forming region ON2N with 
VERA at 14 epochs spanned about two years. The day 
of year (DOY) of these epochs are 112, 206, 246, 305, 350 



in 2006, 009, 051, 095, 132, 222, 277 in 2007, 005, 104, 
and 192 in 2008, which are 53847, 53941, 53981, 54040, 
54085, 54109, 54151, 54195, 54232, 54322, 54377, 54470, 
54569, and 54657 in the Modified Julian Day (MJD), re- 
spectively. The data of the first to third epochs were 
used to estimate the coodinates of H 2 masers. We 
analyze the data of the rest 11 epochs in this study. 
At each epoch, the H2O masers in ON2N and a posi- 
tion reference continuum source, ICRF J201528. 7+371059 
(hereafter J2015+3710), were simultaneously observed 
in a dual-beam mode for about 10 hours. The typ- 
ical on-source integration time was 6 hours for both 
ON2N and J2015+3710. J2015+3710 is listed in the 
VLBA Calibrator Survey 2 (VCS2: Fomalont et al. 2003). 
J2015+3710 was detected with a peak flux density of 0.8- 
2 Jy in each epochs. The separation angle between ON2N 
and J2015+3710 is l.°27. The instrumental phase differ- 
ence between the two beams was measured continuously 
during the observations by injecting artifical noise sources 
into both beams (Honma et al. 2008a). Left-hand circular 
polarization signals were sampled with 2-bit quantization, 
and filtered with the VERA digital filter unit (Iguchi et 
al. 2005). The data were recorded onto magnetic tapes 
at a rate of 1024 Mbps, providing a total bandwidth of 
256 MHz, which consists of 16 IF channels with a band 
width of 16 MHz each. One IF channel was assigned to 
ON2N, and the other 15 IF channels were assigned to 
J2015+3710, respectively. Correlation processing was car- 
ried out on the Mitaka FX correlator. The frequency and 
velocity spacings for ON2N were 15.625 kHz and 0.21 km 
s _1 , respectively. 

Data reduction was conducted using the NRAO 
Astronomical Image Processing System (AIPS). An am- 
plitude calibration was performed using the system 
noise temperatures during the observations. For phase- 
referencing, a fringe fitting was made using the AIPS task 
FRING on J2015+3710 with a typical integration time of 
1 min and a time interval of 30 sec. The solutions of the 
fringe phases, group delays, and delay rates were applied 
to ON2N in order to calibrate the visibility data. Phase 
and amplitude solutions obtained from self-calibration of 
J2015+3710 were also applied to ON2N. Visibility phase 
errors caused by the Earth's atmosphere were calibrated 
based on GPS measurements of the atmospheric zenith de- 
lay which occurs due to tropospheric water vapor (Honma 
et al. 2008b) . After the calibration, we made spectral-line 
image cubes each of which extent on the sky is 1024 x 1024 
pixels with 0.05 mas using the AIPS task IMAGR. The 
typical size of the synthesized beam was 1.2 x 0.9 mas with 
position angle of —50°. The rms noise for each image was 
approximately 0.1-1 Jy beam -1 . The signal-to-noize ratio 
of 7 was adopted as the detection criterion. 

3. Results and Discussion 

3.1. Overall Properties of H 2 Masers in ON2N 

Figure 2 shows the scalar-averaged cross-power spec- 
tra of H 2 masers in ON2N observed with the VERA 
Mizusawa-Iriki baseline on 2007/051 and 2008/104. The 
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intense emissions with the flux density of > 100 Jy were 
detected at the LSR velocity range from — 5 to 5 km s _1 . 
The center of the velocity range of the intense emissions is 
close to the LSR velocity of the associated molecular cloud 
at «lsr = ± 1 km s" 1 of the NH 3 and CS lines (Olmi & 
Cesaroni 1999; Codella et al. 2010). The majority of Wlsr 
of H2O masers is found in the range — 12 < vlsr < 9 km 
s _1 seen in the previous monitoring observations in 1995- 
2004 (Lekht et al. 2006). We also detected the blueshifted 
components at vlsr = — 33 and —17 km s _1 and the red- 
shifted component at 19 km s" 1 , which were not detected 
in the preveous observations (e.g. Lekht et al. 2006). 

Thirty H2O maser spots were detected over the half of 
year and at more than three epochs. They were in the LSR 
velocity range from —33 to 19 km s , and distrubuted 
with an area of 1.0" xO. 4". Figure 3 shows the distri- 
bution of internal motion of the maser spots in ON2N. 
The reference position of the map is set at the position of 
a maser spot at i>lsr = 0.1 km s -1 , which is estimated 
to be (a, <5)j2ooo.o = (20 h 21 m 44?01225, 37°36'37"4844). 
Although the H2O masers are located at 2" south from 
the peak of the 6 cm radio continuum emission (Wood & 
Churchwell 1989), they are spatially coincident with the 
peaks of the 7 mm radio continuum emission and the NH3 
(3,3) emission (Carral et al. 1997; Codella et al. 2010). 

3.2. Parallax and Proper Motion 

The absolute motion of the respective H 2 maser spot, 
i.e. its motion with respective to the position reference 
source J2015+3710, is given by the sum of the proper mo- 
tion and the annual parallactic motion. In order to sepa- 
rate these two motions we performed monitoring observa- 
tions of the H2O masers of ON2N for about two years. 

We made a combined parallax fit, which means a fitting 
of the positions of 30 H2O maser spots to a common par- 
allax, but different proper motion and position offset for 
each spot. Figure 4 and Table 1 show the results of the 
combined parallax fit. As can be clearly seen in Figure 4, 
the observed points demonstrate a sinusoidal modulation 
with a period of 1 year caused by the annual parallax. 
For this fitting, we assigned independent "error floors" in 
quadrature with the formal position fitting uncertainties. 
Trial combined fits were conducted and a separate reduced 
X 2 (per degree of freedom) statistics was applied for the 
right ascension and declination residuals. The error floors 
of 0.088 mas and 0.111 mas in the right ascension and the 
declination respectively, were then adjusted iteratively so 
as to achive a reduced \ 2 P cr dgrcc of freedom near unity 
in each coordinate. Combining all of the fittings, we ob- 
tained the trigonometric parallax of the H2O maser spots 
0.261 ±0.009 mas. The parallax gives the heliocentric dis- 
tance of ON2N 3.83 ± 0.13 kpc. 

In Table 1, we also show the estimated parallaxes using 
individual fitting for each maser spot. We made this indi- 
vidual fitting only for 14 maser spots which were detected 
over one year. The obtained parallaxes from the individ- 
ual fitting arc consistent with each other and with the 
result of the combined fit. This means that the parallax 
obtained by the combined fitting is reliable. 



The systemic motion of the source can be estimated as 
an averaged motion of all maser spots, provided that in- 
ternal motion is random or symmetric. We believe that 
this may be reasonable for ON2N because of the follow- 
ing two reasons. The averaged radial velocity of all maser 
spots is —1.2 km s _1 which is close to the systemic radial 
velocity derived from the associated molecular cloud. This 
suggests that the maser spots move rather symmetrically. 
Figure 3 shows the residual proper motion vectors, which 
are differences between the individual proper motions and 
the average. We did not find any strong asymmetric mo- 
tion. Therefore, we find that the reasonable absolute 
proper motions are not biased by the internal motions 
of ON2N. Thus, the systemic proper motion of ON2N is 
estimated to be (fj, a cos5,fi s ) = (-2.79±0.13,-4.66±0.17) 
mas yr _1 as the average motion of all 30 maser spots. 

Using the Galactic coordinates of ON2N (l,b) = 
(75.°78, — 0.°34), the proper motion components in the 
galactic coordinates are calculated to be (fxi cosb, fib) = 
(—5.42 ± 0.16, —0.36 ± 0.14) mas yr _1 corresponding to 
the linear velocity of (v u v h ) = (-98.4 ±2.9, -6.6 ±2.6) km 
s _1 at the distance of 3.83 kpc. We note that these values 
are still affected by the solar motion, because the observed 
proper motion is not relative to the LSR but to the Sun. 

To convert this observed velocity to that with re- 
spect to the LSR, we have to fix the solar motion rel- 
ative to the LSR. As mentioned in §1 we use the solar 
motion in the traditional definition of (Uq,Vq,Wq) = 
(+10.3, ±15.3, +7.7) km s" 1 after Reid et al. (2009a) 1 . 
We note that this set of values is very close to the one 
found by Miyamoto & Zhu (1998) (see Table 5) on the 
basis of the stellar motion analysis of the HIPPARCOS 
proper motions in the solar neighbourhood. Using this 
traditional solar motion, the observed proper motion is 
converted to (^ cos 6, /i b ) = ( - 5 . 76 ± . 1 6 , . 02 ± . 1 4) mas 
yr _1 and the corresponding linear velocity is (17, v/,) = 
(-104.6 ±2.9, 1.1 ±2.6) km s" 1 relative to the LSR. 

3.3. Derivation of the Galactic Constants 

Based on the observed LSR velocity of ON2N we believe 
that the source is located at or close to the solar circle. 
The small proper motion along the galactic latitude, Vb = 
1.1 ± 2.6 km s _1 , shown in §3 supports that ON2N rotates 
circularly around the Galactic center. 

For an object on the solar circle, the galactocentric dis- 
tance of the Sun or that of the source, Rq, is estimated 
from the heliocentric distance of the source as 

Ro = -^-j, (1) 
2 cost 

where I is the galactic longitude of the source (see Figure 
1). Our estimation of the heliocentric distance of ON2N 
yields D = 3.83 ± 0.13 kpc. It gives Rq = 7.80 ± 0.26 kpc, 
if ON2N is exactly located on the solar circle. This value 
is close to the previous estimations; Rq is estimated to be 

1 This velocity is consistent with the provisional solar motion we 
used, but the 3 dimensional linear velocity (Uq,Vq,Wq) = 
(+10.0, +15.4, +7.8) km s _1 shown in Kerr k. Lynden-Bell 
(1986) is inconsistent to the definition by themselves. 
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Table 1. The obtained values of parallax n and proper motions ^ QCOS i5 and fig for H2O maser features in ON2N. 
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Combined fit 0.261 ± 0.009 

Average -2.79±0.13 -4.66±0.17 



Column (3), (4): Right ascension and declination offsets relative to the positon of the maser spot at iJLSR = 0.1 km 
s -1 , and (a,5)j2ooo.o = (20 h 21 m 44?01225,37°36'37f / 4844). 

Column (5): Each alphabetical letter represents the epoch with maser detection. A, B, C, and K mean the 11 
epochs from 53847, 53941, and 54657 in MJD, respectively. A dot represents the epoch without the detection. 
Column (6): Parallax estimated from the individual fitting. 

Column (7), (8): Motions on the sky in the directions along the right ascension and declination. 
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Fig. 2. A scalar-averaged cross-power spectra of H2O masers in ON2N observed with the VERA Mizusawa-Iriki baseline at 2007/051 
and 2008/104. 



8.0±0.5 kpc from a combination of many methods reviewd ^£ 
by Reid (1993), 7.9±g;| kpc fr om a parallax measurement Ro 
of H 2 masers in Sgr B2 with VLB A (Reid ct al. 2009b), 

and 8.28 ±0.44 kpc using the orbits of stars around Sgr = — a °/ i; > v r 

A* from VLT and Keck data (Gillessen et al. 2009). 

For a source on the solar circle with purely circular ro- 
tation, its proper motion velocity along the galactic lon- 
gitude, vi , gives the galactic rotation velocity of the LSR 
or that of the object as 
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(5) 



e 
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where I is the galactic longitude of the source. Our mea- 
sured value of vi = —104.6 ±2.6 km s _1 gives 9o = 213 ±5 
km s _1 . This value is smaller than that estimated by Reid 
et al. (2009a), 6 = 254 ± 16 km s -1 but close to the IAU 
recommended value of 9o = 220 km s _1 . 

Our derivations of the Galactic constants are strongly 
dependent on the assumption of the location of ON2N 
in the MWG. However, we found that the ratio of the 
Galactic constants, ©0/-R0, which is the anglular velocity 
of the LSR, tto, can be estimated less dependently on the 
assumption. The ratio O0/-R0 can be estimated, even if 
ON2N is not exactly located on the solar circle, but near 
there. 

For a source on pure circular rotation at any position in 
the galactic disk, its radial and tangential velocities with 
respect to the LSR can be written as 



where do is a constant to convert the unit from an angular 
velocity to a linear velocity and ao = 4.74 km s _1 mas -1 
yr kpc -1 . For a source near the solar circle, its v r is nearly 
zero. In this case, equation (5) yields 

ctoW, (6) 



(2) ^ 



v, 



Oo 
Rq 
Oo 



Rq sin/, 



e 



Rncosl — —D. 
R 
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where R is the actual galactocentric distance of the source, 
and D is the heliocentric distance of the source, and O is 
the galactic rotation velocity of the source. Equations (3) 
and (4) yield 



#0 

which is free from Rq. Actually we found that Qq/Rq 
is a nearly constant at 6 < Rq < 10 kpc and the ratio is 
obtained as ©0/-R0 = 27.3 ±0.8 km s -1 kpc -1 using the 
D = 3.83 ±0.13 kpc, m = -5.76 ±0.16 mas yr" 1 , and 
v r = 0± 1 km s _1 . This value is close to the value of 
Q0/R0 = 28.7±1.3 km s _1 kpc -1 obtained from a tangent 
point source, ONI (Nagayama et al. 2011), and the value 
of Oo/i?o = 28.6 ±0.2 km s _1 kpc -1 obtained from the 
proper motion measurement of Sgr A* (Reid & Brunthaler 
2004), which is revised using the traditional definition of 
the solar motion by us. However, this value is inconsistent 
to that derived from the IAU recommended values 220 
km s _1 / 8.5 kpc=25.9 km s^ 1 kpc -1 . This estimation 
gives another constraint on the Galactic constants which 
is independent of the Oort constants derived from stellar 
motion near the Sun. 
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Fig. 3. (a): Distributions of H2O mascrs (color filled circles) superimposed on the 6 cm radio continuum emission (cyan contour; 
Wood & Churchwell 1989), the 7 mm radio continuum emission (gray contour; Carral et al. 1997), and the NH3 (3,3) emission (black 
contour; Codclla et al. 2010). (b): Internal motion vectors of H2O masers. The spot color shows the LSR velocity. The arrow at 
the top left corner shows a internal motion of 1 mas yr _1 , corresponding to 18.2 km s — 1 at a distance of 3.83 kpc. (c): Close-up to 
the central part of (b). 
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Fig. 4. Parallax of the H2O mascrs in ON2N. The data for the different maser positions are slightly shifted in time for clarity. 
Individual proper motions and position offsets are removed. The left and right panels show the annual parallax in right ascension 
and declination, respectively. The numbers in each panel show the ID number of the spot listed in Table 1. 
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